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We investigate asymptotical behavior of numbers of long Hamiltonian walks (HWs), i.e. self- 
avoiding random walks that visit every site of a lattice, on various fractal lattices. By applying 
an exact recursive technique we obtain scaling forms for open HWs on 3-simplex lattice, Sierpinski 
gasket, and their generalizations: Given-Mandelbrot (GM), modified Sierpinski gasket (MSG) and 
n-simplex fractal families. For GM, MSG and n-simplex lattices with odd values of n, the number 
of open HWs Zn, for the lattice with N 1 sites, varies as w N N J . We explicitly calculate the 
exponent 7 for several members of GM and MSG families, as well as for n-simplices with n = 3,5, 

and 7. For n-simplex fractals with even n we find different scaling form: Zn ~ oJ n , where df 
is the fractal dimension of the lattice, which also differs from the formula expected for homogeneous 
lattices. We discuss possible implications of our results on studies of real compact polymers. 
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I. INTRODUCTION 

Self-avoiding walks (SAWs) have long been used in the 
studies of configurational statistics of polymer chains in 
solution Due to excluded volume effect, at high tem- 
peratures T (good solvent) long polymer chains are in 
swollen configurations. At low temperatures (poor sol- 
vent) polymers are in a collapsed state, caused by the 
attractive interactions of different sections of a polymer, 
mediated by a solvent. The transition between these two 
states occurs at the so-called 9 temperature, at which 
excluded volume and attractive forces balance. Whereas 
the swollen and phases has been well investigated by 
now, the entropic scaling of the collapsed phase is still 
an open issue. A closely related problem is the scaling 
of Hamiltonian walks (HWs), which are SAWs that visit 
all the sites of the underlying lattice. HWs are believed 
to represent the T = limit of collapsed polymers, and 
they are also used in the studies of polymer melting 0, 
as well as in the context of protein folding [3| . 

The number Zn of HWs on homogeneous lattices with 
N S> 1 sites is expected to behave as 



connectivity constant, defined as 



u N vfN a . 



(1.1) 



Here a — (d — 1) / d, where d is the dimensionality of the 
lattice, fis is some constant less than 1, and to is the 
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Proposed scaling form for HWs differs from the ordinary 
SAW case (swollen polymer), where average number of 
iV-step SAWs, for large N, behaves as uj n N a , and where 
the critical exponent a depends only on d (which is not 
the case for HWs). The term fj,g in the HW case is 
expected on the basis of the exact study of HWs on the 
Manhattan lattice Q, as well as on the conjecture that 
collapsed polymer (globule) has a sharp boundary, so 
that a surface tension term should arise la] . The scaling 
form (jl.ip was confirmed by Owczarek [fj for collapsed 
partially directed SAWs on the square lattice. Baiesi et 
al p} recently performed extensive Monte Carlo simula- 
tions, which gave strong evidence that (| 1 . 1 [) is also sat- 
isfied for undirected collapsed SAWs on the square lat- 
tice. There are hardly any results for higher dimensional 
lattices. To the best of our knowledge, the only clear 
indication that for collapsed SAWs on three-dimensional 
lattices there exists a surface term, as predicted by p. II) , 
was obtained by Grassberger and Hegger Q via Monte 
Carlo simulations. There are also no results concerning 
scaling form for the collapsed SAWs on disordered lat- 
tices. Having all that in mind it might be useful to study 
HWs on fractals. 

Fractal lattices are somehow intermediate between ho- 
mogeneous and disordered ones, and their hierarchical 
and scale invariant structure often allows an exact recur- 
sive treatment of various physical phenomena. The or- 
dinary SAW model (corresponding to the polymer chain 
at high temperatures) has been studied extensively in 
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the past on different fractals These studies con- 

tribute to a better understanding of how dimension and 
topological structure of underlying space affect the criti- 
cal behavior of SAWs in general. There are fewer papers 
about self- interacting SAWs on fractals [2l[- 25 1, and nei- 
ther of them establishes the explicit form of the partition 
function at low temperatures. The closely related prob- 
lem of finding the number of walks in the limiting HW 
case was analyzed only for closed walks on some fractals 
(26l [27| . Here for the first time we investigate the scaling 
forms for the number of open HWs on a set of fractal 
lattices. 

The paper is organized as follows. In Sec. [TT] we de- 
scribe the technique used for an exact enumeration of 
closed HWs on the cases of 3-simplex lattice and Sierpin- 
ski gasket fractal, and extend it for enumeration of open 
HWs on these lattices. In Sections Mil and HVl this tech- 
nique is generalized for enumeration of HWs on Given- 
Mandelbrot (GM) and modified Sierpinski gasket (MSG) 
families of fractals. We find that there is no surface term 
in the scaling form for the number of HWs on neither of 
these lattices. For the number of closed HWs on MSG 
fractals we show that it varies as u> N , which is the same as 
in the case of GM fractals [27| . For open HWs we obtain 
scaling form uj n N" ( for both fractal families, with ui and 
7 depending on the particular parameters of the fractals. 
Furthermore, exponent 7 does not have the same value 
for corresponding fractals from these two families, as is 
the case for ordinary SAWs [l(J HH . Similar violation of 
the universality was found for two-dimensional (2d) ho- 
mogeneous lattices [1^] . We extend our analysis to higher 
dimensional n-simplex lattices in Sec. |Vl where general 
exact scheme for obtaining numbers of HWs is described. 
It turns out that scaling form for HWs strongly depends 
on the parity of the fractal parameter n. In particular, 
for lattices with odd values of n we find the same scaling 
form as for HWs on GM and MSG fractals, whereas for 
even n we get Zn ~ oj n fi N ' 1 , with df being the frac- 
tal dimension of the lattice, and different values of \i for 
closed and open HWs. All obtained results are summa- 
rized in Sec. VI, and some technical details are given in 
the appendixes. 



II. HAMILTONIAN WALKS ON 3-SIMPLEX 
LATTICE AND SIERPINSKI GASKET FRACTAL 

In order to explain the method used in this paper for 
obtaining the scaling forms of numbers of both closed 
and open HWs we start with the simplest case: HWs 
on the 3-simplex lattice, and proceed with HWs on the 
Sierpinski gasket. 



A. 3-simplex lattice 

To obtain the 3-simplex lattice 30] one starts with a 
complete graph of three points and replaces each of these 




(a) 



(b) 




FIG. 1: (a) An open HW on a 3-simplex of order I = 4. The 
first (b) and the second step (c) of the coarse-graining process. 
Gray triangles in (b) and (c) represent 3-simplices of order 
I — 2 and I — 3, respectively, whereas curved lines correspond 
to the coarse-grained walk. Different types of configurations 
within the third and the second order simplices are encircled, 
(d) Open HW on 3-simplex of order (I + 1), consisting of two 
A-type parts and one -B-type part within the corresponding 
simplices of order I (gray triangles). 



points by a new complete graph of three points. The sub- 
sequent stages are constructed self-similarly, by repeat- 
ing this procedure. After I such iterations one obtains 
3-simplex of order I (see Fig. da)), whereas the com- 
plete 3-simplex lattice is obtained in the limit I —> 00. 
The number of vertices within the Zth order 3-simplex 
is AT; = 3', and the fractal dimension of the lattice is 
d f = In3/ln2. 

In Fig. [UJa) we give an example of open HW on a 3- 
simplex of order I = 4. Performing a coarse-graining 
process one notices in Fig. [IJb) that this walk can be de- 
composed into parts corresponding to the second order 
simplices, which have two possible configurations: one 
traversing the simplex (B-type configuration) and the 
other consisting of two strands, one of which is traversing 
the simplex, and the other with one end at a corner vertex 
of the same simplex, and the second end anywhere within 
it (C-type configuration) . The coarse-graining can be ap- 
plied once more (Fig. QIc)), leading to a HW consisting 
of three parts corresponding to three simplices of order 
/ = 3, two of them being of the type B, and the third 
consisting of two strands ( "legs" ) , each of them having 
one end-point at the corner vertex of the corresponding 
simplex and the other within it (D-type configuration). 
In a similar way, it is not difficult to see that each open 
HW on the (/ + l)th order simplex can be decomposed 
in two ways, either as in the sample shown in Fig. [IJc), 
or as in Fig. HJd), with one -B-type part and two ^4-type 
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FIG. 2: (a) Four possible configurations of A-type HWs on 
a 3-simplex of order (I + 1), with the corresponding numbers 
of HWs. Gray triangles represent the ith order simplex, (b) 
Two possible A-type HWs on the first order 3-simplex. Small 
circles represent sites the lattice consists of. 



parts (one strand with one end at the corner vertex and 
the other end anywhere within the corresponding Zth or- 
der simplex) . If we denote the number of open HWs on 



the (I + l)th order simplex by Z Q 
relation is valid: 



then the following 



= 3B<'> 



where A {1 \ 5« and D<0 are the numbers of HWs of 
types A, B and D within the Z-th order simplex, respec- 
tively. Decomposition of closed HWs is even simpler, 
leading to the recursion relation: 



AW) 



(2.2) 



where Z. 



(z+i) 
c 



is the overall number of the closed HWs 



on the simplex of order (I + 1). Number C*(0 of C-type 



HWs does not contribute directly to Zq +1 ^ and Z { 
but one can obtain (see Fig. [5]) the following recursion 
relations: 



Al+i) 



A (t+V 

£)(* + !) 



2B 2 (A + C), 
B\ 

B 2 (A + 3C), 
B{AAC + 2A 2 + 3BD 



(2.3) 
(2.4) 
(2.5) 

6C 2 ) , (2.6) 



for the numbers of HWs corresponding to simplices of 
two successive orders (where we have omitted index (I) on 
the right-hand side of the relations). Iterating these rela- 
tions, starting with the initial values A^ 1 ' = 2, B^ 1 ' = 1, 
CW = 1, and I)W = 2, one can calculate A^, B®, 
C«, and for any I, and, consequently, the numbers 



and Z, 



(i+l) 



Furthermore, these relations can be 



exactly analyzed, as follows. From (|2.4p . and the cor- 
responding initial condition, trivially follows B^ = 1, 
which also transforms the set of relations (|2.3|) i (12. 5|) for 



one-leg walks A and C into a simple system of difference 
equations: 

A d+D = 2{A (D + C W) ( C d+D = A {i) + 3C (i) } 
whose solution is 

A w = i(4' + 2), C il) = i(4'-l). (2.7) 



Then, from recursion relation (J2T6J) straightforwardly fol- 
lows the difference equation for the number of two-leg 
walks D: 



£,(1+1) = ± 4 2r + 3D (l) 

3 



whose solution is 



n(0 — — I6 r + — 3 r - - 
^ - 39 ib + 13 d 3 



(2.8) 



Finally, putting (|2~7) . (|2^|) . and B<0 = 1 in formulas for 
the overall numbers of open (|2.1[) and closed (|2.2[) HWs 
one gets 



m 25 ,1,3,1 
= 16T39 16i + 3 4 ' + 13 3 ' + 3 ' 



7 (0 



1 



(2.1) From the last equation trivially follows that connectivity 



constant is u> = 1 and Zq = u> Nl for every I, whereas the 
number of open HWs for I ;§> 1 behaves according to the 
asymptotic formula 



Z, 



(i) 



(2.9) 



where Ni — 3 l is the number of vertices of the 3-simplex 
of order I, and 7 = Inl6/ln3 = 2.52372 



B. Sierpinski gasket 

In a similar way, one can analyze HWs on the Sierpin- 
ski gasket (SG) lattice. SG is a well known fractal lattice, 
which can be constructed recursively, starting with the 
generator (gasket of order r = 1), which consists of three 
unit equilateral triangles, arranged to form a twice larger 
triangle (see Fig. OJa)). The subsequent fractal stages 
are constructed self-similarly, by replacing each of the 
unit triangles of the initial generator with a new gener- 
ator. To obtain the rth-stage fractal lattice (rth order 
gasket), this process of construction has to be repeated 
(r — 1) times, and the complete fractal is obtained in 
the limit r —> 00. The numbers of sites on the rth or- 
der gasket is equal to N r = |(3 r + 1). SG resembles 
3-simplex lattice and indeed has the same fractal dimen- 
sion df = In 3/ In 2. An open HW on a third order gas- 
ket is shown in Fig. [3l together with its coarse-grained 
versions. Comparing with Fig. [1] one can observe that 
larger number of types of possible HW configurations ex- 
ists on SG than in the case of 3-simplex lattice. There 
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FIG. 3: (a) An open Hamiltonian walk on Sierpinski gas- 
ket of order r = 3. Vertices corresponding to the end-points 
of the walk are colored black, for the sake of easier recogni- 
tion. The first (b) and the second step (c) of coarse-graining 
process. Gray triangles in (b) and (c) represent gaskets of 
order r = 1 and r — 2, respectively, whereas curved heavy 
lines correspond to the coarse-grained walk. Different kinds 
of configurations within the first and the second order gaskets 
are encircled. 



are exactly 8 different types of walks, and they are de- 
picted in Fig. 0] The corresponding numbers of HWs on 



the rth-order gasket will be denoted by A^ r 
B«, B<f\ C«, DW, and D{ r \ 
the following recursion relations: 



Ar) 



A. 



(r) 



These numbers fulfill 



A {r+1 '> = 2(B 2 C + ABB 1 +A 1 B 2 



(2.10) 



A { { +1) = AB^ + 2A 1 BB 1 +2A 2 B 2 + 2BB X C ,(2A1) 



4 +1) 

£(r-+l) 
C (r+1) 



A\Bi + 2A 2 BBi + Bi C , 



2B 2 B 1 



B[ r+1) = 2BB 



D 



(r+l) 



(2.12) 
(2.13) 

AB 2 + 2B(AiBi + A 2 B + 3£?iC), (2.14) 
AAB X C + 4ABA 2 + ±A Y BC + ABBiD 
AAB l A l +2B{ZC 2 + 2BD±+ A\), (2.15) 

AA Y A 2 B + 2AA 2 B Y + 2A 1 CB 1 + ABB 1 D 1 
AA 2 CB + ZC 2 B 1 +A 1 2 B 1 +B 1 2 D, (2.16) 




FIG. 4: Possible types of HWs on rth order gaskets (gray 
triangles), needed to decompose any HW on gasket of order 
r + l. 



However, this is a tedious task, since all of these numbers 
quickly become extremely large. To avoid that, we ana- 
lyze recursion relations in a similar manner as was done 
for 3-simplex lattice. 

Introducing new variables: x r = A^/B^ r \ x\ r = 



Vr 



AF/B<r\ x 2r = A\ r) /B^ 
DW/flW, z lr = dP/bM 



= C^/B^, z r = 

and noticing that 

Bi'/B^ — 2/3 is exactly satisfied for any r, from 
(|2.10j) - (|2.16j) one gets the new recursion relations: 

/ x' 




/I 
l 

3 



\\ 

2 

_L 

3 



1 
_L 

3 

3/ 



( x\ 

Xi 

x 2 
\ V 1 

h 
h 



(2.18) 



(2.19) 



where 
fx 
h 



3 2 9 2 

2xx\ + —x 1 + 3xx 2 + 2xy + 3xiy + —y , 

\^x\ + xx 2 + 3x!X 2 + x x y + 3x 2 y + | y 2 . 



From (|2 . 18[) for r 3> 1 follows x r , x\ r , x 2rl y r <~ A r , 
where A = 4 is the only eigenvalue of the 4x4 matrix of 
this recursion relation which is larger than 1, and then, 
from (rZHD, we find z r , z lr ~ 4 2r . 

To obtain the asymptotic formula for the number of 
open HWs, we express Zq from (|2.17p as 



8 / n \3 

= 3V J ^ XlrX2r + Zl r) ■ 



with the initial values: = 4, a{ 1] = 6, A 2 1] = 3, 

SW = 3, fif 3 = 2, C« = 4, DW = 6, and = 5. 
Iterating recursion relations (|2.10|) - (|2.16p one can calcu- 
late numbers of all possible types of HWs, in principle 
for any r, and then eventually find the overall number 



Since in [27} it was exactly shown that B^ = const uj Nr , 
with uj — 12 1 / 9 , one finally gets the scaling form: Zq^ ~ 
ui N rN?, with 7 = Inl6/ln3 = 2.52372.... For the 

(r) 

sake of comparison, we have generated Z Q by direct 
numerical iteration, for r up to 100. In Fig. [5]we depict 



of open HWs on the (r + l)th order gasket, by {Y\\Z y ')/N ri as a function of (lnN r )/N r , for 21 < r < 30 



putting them into the formula: 
Z ( o +1) = 12 A^ r) A 2 {r] Bx (r) + 6 (Bx (r) ) 2 L>x (r) . (2.17) 



As one can see, value of ()d.Zq)/N t for this range of r is 
quite close to In a;, whereas by fitting these data to linear 
function one obtains approximate value: 7/,-j = 2.3. For 
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FIG. 5: Values of (\n Z Q ) / N r , found by direct numerical 
iteration of relations pl0]l - (|2lg]l and pTTTjl . for 21 < r < 30. 
Line connecting the points has been obtained by linear fitting, 
and uj = 12 1/9 . 



ing the smallest triangles with the generator, and then 
repeating this procedure recursively ad infinitum, one ob- 
tains fractal lattice characterized with the integer 6. For 
b = 2, 3, ...,oo the complete so called Given-Mandelbrot 
(GM) family of fractals is obtained [3l[. SG is the first 
member of this family, with the scaling parameter 6 = 2. 
The fractal dimension of GM fractal with scaling param- 
eter b is equal to df = ln[6(6+l)/2]/ ln6, and the number 
of sites N r at the rth stage of fractal construction is 



N r 



6 + 2 



1) 



1 



6 + 2 



(3.1) 



The overall number of open HWs on the (r + l)th stage 
of construction of any GM fractal can be expressed in 
terms of numbers of 8 HW types within the rth order 
stage, in a similar manner as in the case of SG (Fig. [1]). 
Recursion relations for numbers of B, and £?i-type walks 
on two successive stages of fractal construction have the 
following form: 



B' = pB- 



1 B 1 



b-l 



010— 1) i 

B[=pB^B 1 \ (3.2) 



larger number of iterations "ffu becomes closer to the ex- 
act value of 7. For instance, fitting the numerical data 
for 91 < r < 100 gives (7 - 7/«)/7 = 2%. 

Scaling form obtained for the number of open HWs on 
SG is the same as the one found for the 3-simplex lat- 
tice. As was shown in [27j, the numbers of closed HWs 
on SG scale as ui Nr , again the same as in the case of 3- 
simplex. Equality of the exponents for SG and 3-simplex 
is in accord with the fact that SAWs on these two lattices 
belong to the same universality class 28]. On the other 
hand, it is known that exponents for HWs on different 2d 
Euclidean lattices have different values [2i|, which is ex- 
plained to be a consequence of the frustration, induced 
by the strong constraint that all the sites must be vis- 
ited. It is believed that a relevant physical measure of 
this frustration is the number of contacts per monomer, 
i.e. vertex pairs which are not adjacent along the HW, 
but are the nearest neighbors on the lattice. Neverthe- 
less, the number of contacts on 3-simplex is one, whereas 
it is two on SG, so that one could have expected differ- 
ent values of 7. In order to gain a deeper insight into 
the problem of universality and frustration of HWs on 
lattices embedded in 2d space, in the next two sections 
we analyze the asymptotic behavior of HWs on the ap- 
propriate generalizations of 3-simplex and SG fractals. 



III. HAMILTONIAN WALKS ON 
GIVEN-MANDELBROT FRACTALS 



One possible way to generalize the SG fractal is to 
start with a generator that consists of 6(6 + l)/2 unit 
equilateral triangles, arranged to form a 6 times larger 
triangle. Enlarging the generator 6 times and substitut- 



as was shown in [27|, whereas the numbers of one- leg 
configurations: A, A\, A2, and C, satisfy a closed set of 
recursion relations, which can be put in matrix form as: 



(A! 

A 

\ C 



( an au 013 an 

»21 «22 023 G24 

»31 a 32 Q-33 a 34 

\ »41 «42 fl43 «44 




(3.3) 



In Ap- 



where a-y are polynomials in B^ and B^ . 
pendix A we prove that each of these polynomials has 
only one term, which is of the form -[B^] 13 ^^] 6 , 
where (3 + 6 = 6(6 + l)/2 - 1. The coefficients K tj de- 
pend only on 6, and each of them can be expressed in 
terms of the number p (appearing in equations (|3.2p ). 



and additional 8 numbers: a, a%, a2, c, a', a' l7 a 2 , d '. Intro- 



ducing new variables 



j- 2 



<7BW, y r = CW/B W , relation 



gives the 



new recursive relation: 
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"1 


a-2 
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tp 


t 2 p 


tp 


t(2a-p) 


<n 


»2+P 


c 




2p 
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tp 
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i 2 (a-p) 


ai t 


a 2 +2p 


c I 




2p 


2p 


2p 


2p 




a' 1 


'A 


"2 
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tp 


p 





( x \ 




Xi 








\ y J 



V 

(3.4) 

where t = B^/B^ is a coefficient that depends only 
on 6, as can be seen from (13. 2j) . The particular values 
of the coefficients appearing in the last 4x4 matrix, for 
2 < 6 < 7, are given in Appendix A. For all 6 considered 
here, this matrix has only one eigenvalue A larger than 
1, implying that for r>l parameters x r , x\ r , X2 r , and 
y r grow as A r , and consequently, numbers A^ r \ A\- r \ 
^2^, and asymptotically behave as X r B^ r \ 
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Numbers of two-leg HWs satisfy recursion relation: 

(&)-(£ £)U) + U)- ( "> 

where bij are polynomials in and of the power 
b(b + l)/2 — 1 (see Appendix B), whereas F and Fi 
are of the form: cIaaA 2 + cIaa 1 AA 1 + dAA 2 AA 2 + 
dAcAC+dA 1 A 1 A 2 1 +dA 1 A 2 A 1 A 2 + d Al cA 1 C+dA 2 A 2 Al + 
dA 2 cAiC + dccC 2 , with dxy being polynomials in B^ 
and B± of the power b(b + l)/2 — 2. With new param- 
eters defined as z r = D^/B^ and z lr = D^/B^, 
from (|3.5p follows simpler recursion relation: 

f"' Wl? U mi2 ]f Z ) + f {), (3.6) 
y zi y y m 2 l TO 22 J \Zl J \h J 



where rriij are numbers, which do not depend on r, but 
only on b, whereas / and /i are polynomials in x r , X\ r , 
x 2r and y r . In particular, it can be shown that tou = 
to + 1, mi2 = 2n/t, m 2 \ — mt/2, m 22 = n + 1, with 
rn = fe(fe — l)/2, and n = b — 1, whereas / and /1 have 
the form: a^x 2 + a XXl xxi + a XX2 xx 2 + a xy xy + a XlXl x\ + 
xiy + a X2X2 x\ + a X2V x 2 y + a yy y 2 , where 
dxy are again numbers depending only on the scaling 
parameter b. It turns out that for all b considered here A 2 
is larger than eigenvalues A^> of the matrix to^ , implying 
that z r , z\ r ~ A 2r . 



The number of open HWs on the gasket of order 

(r + 1), can be expressed as: 



Z% +1) = kiAAiB m ~ 3 B™ +2 + k 2 AA 2 B m - 2 B™ +1 + ksA^B^B? + k 4 ACB m - 3 B™ +2 
+ k 5 A 1 CB m - 2 B1 +1 + k 6 A 2 CB m - l B^ + k 7 A 2 B m - i B™ +3 + k 8 A 2 B rn - 2 B^ +1 
+ k 9 A 2 2 B m B 1 ^- 1 + k w C 2 B m - 2 B? +1 + kuDB m - 2 B™+ 2 + k^DiB" 1 ^ 1 B" +1 , (3.7) 



where we have suppressed index r on the right-hand side 
of this relation, and ki are numbers that depend only 
on b. Substituting established asymptotical behavior of 
Ai, C and Di in the latter expression, one finds that all 
terms on the right-hand side of equation (|3.7p have the 
same asymptotical form, so that 



Z. 



(r+l) 



A 2r (B w ) ! 



Since in [27j it was shown that numbers B^ and B[ r ^ 
for large r behave as cu Nr , where the number of sites N r 
is given by (|3.1|) . it follows that 



Z. 



(r) 
O 



Lo" r N? , with 7 



lnA(fr) 

faQ+i) 
2 



hi 



(3.8) 



The values of 7 for 2 < b < 7 are equal to 
2.5237..., 2.1841..., 2.3411..., 2.2461..., 2.2981..., 
and 2.2755 . . ., respectively. One should mention here 
that number of closed HWs asymptotically behaves as 
Z^ ~ to Nr , as was established in (2?| , where also a closed 
formula for lu was derived. 




FIG. 6: (a) The (r + l)th order generator of the GM frac- 
tal with b — 4. Grey triangles represent its 10 constitutive 
elements, each of them being the rth stage of construction 
of the fractal (rth order generator). Most of the corner ver- 
tices (black circles) of the rth order generators belong to more 
than one triangle, (b) The (r + l)th order generator of the 
MSG fractal with b — 4 is obtained by connecting rth order 
generators (grey triangles) via infinitesimal junctions. Curved 
lines represent two similar £?-type HWs which in (a) case con- 
sists of B and B\ type strands within the rth order triangles, 
whereas in (b) case each strand within the rth order triangle 
is of B type, implying that all vertices within the triangle are 
visited (including the third corner vertex). 



IV. MODIFIED SIERPINSKI GASKET 
FRACTALS 

Intrigued by the fact that scaling relations for the num- 
bers of HWs, as well as the corresponding values of the 
critical exponents, are the same for 3-simplex and SG lat- 
tices (Sec. II), here we examine HWs on modified Sierpin- 
ski gasket (MSG) fractals. In particular, the only differ- 
ence between GM fractal and corresponding MSG fractal 



(with the same value of scaling parameter b) is that the 
generator of the MSG fractal is obtained by arranging 
b(b + l)/2 unit triangles in such a way that vertices of 
neighboring triangles are connected with an infinitesimal 
junction, i.e. not glued, as in the GM case (see Fig. [6]). 
This insertion of junctions simplifies recursive scheme for 
counting HWs: instead of eight different types of walks 
in the case of GM fractals, one should consider only four 
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TABLE I: Coefficients q, a, c, a' , and c', appearing in recur- 
sion relations (|4.1|l . and (|4,4[l . for the numbers of B, A, and 
C-type HWs, found by direct computer enumeration of all 
possible corresponding configurations, together with the val- 
ues of connectivity constant ui and exponent 7 for open HWs 
on MSG fractals with 2 < b < 8. 



6 


<7 


a 


c 


a' 


c' 


Id 


7 


2 


1 


2 


2 


1 


3 


1 


2.52372 


3 


2 


6 


10 


4 


8 


1.04729 


2.12343 


4 


4 


24 


42 


20 


40 


1.05269 


2.3816 


5 


16 


134 


228 


116 


196 


1.06824 


2.23644 


6 


68 


932 


1460 


886 


1408 


1.07286 


2.32218 


7 


464 


8656 


12524 


8372 


11604 


1.07875 


2.27302 


8 


3838 


101612 


133764 


103258 


133428 


1.08177 


2.29907 



types of walks in this case. 

Any closed HW within the (r+ l)th order generator of 
MSG can be decomposed into 6(6 + l)/2 strands, each of 
them being a -B-type HW within one of the 6(6 + 1)/2 rth 
order triangle inside it. Here, S-type HW has the same 
meaning as in the case of 3-simplex lattice, i.e. it is a HW 
that enters a generator at one of its corner vertices, and 
leaves it through another one, meanwhile visiting all the 
remaining vertices, including the third corner vertex (see 
Fig. [6J. Recursion relation for the numbers of S-type 
walks is of the form 



B (r+1) = q f B (r)y 



(4.1) 



where q depends only on 6, and, consequently, for r >> 1, 
the number behaves as: 



B {r) 



(4.2) 



The number of closed HWs on the (r + l)th order gener- 
ator Zq is equal to: 



KH-l) 



= const ( B (r) 



meaning that has the same asymptotic form as B^ r \ 
Since the number N r of sites on the rth order generator 
is equal to N r = 3[6(6 + l)/2] r , for large r one gets the 
following scaling form: 



Z% ] ~ B {r) ~ u Nr , with to = qW^m^J . (4.3) 

Values of q and ui for 2 < b < 8 are given in Table |U 

The number Zq +1 ^ of open HWs on the (r + l)th order 
generator is equal to: 



?(r) 



B— 2 (fciA 2 + k 2 AC + k 3 C 2 + k 4 BD) 



where A = A^ r \ B = B {r \, C = and D = have 
the same meaning as in the case of 3-simplex lattice, and 



ki are some constants depending only on b. Numbers A 
and C for two consecutive stages of MSG fractal con- 
struction satisfy the recursion relation: 



(4.4) 



where a, c, a', and c' depend only on 6, and can be found 
by enumeration of the corresponding one-leg HWs (see 
Table|T|. Dividing this relation by recursion relation (14. 1|) 
for the B numbers, and introducing new variables: x r = 
A^/B^ r \ and y r = C^/B (r \ one gets 





(4.5) 



implying that for r> 1 parameters x r and y r behave as 
A r , where 

A = ^- (a + c' + y/(a- c') 2 +4ca') (4.6) 

is the larger eigenvalue of the 2x2 matrix in (|4.5|) . 
The number of two-leg HWs transforms as: 

D 1 = B^- 2 (d AA A 2 + d AC AC + d cc C 2 + d D BD) , 

which, with z r = D^/B^ can be rewritten as 

z r +i = - (d A A x 2 r + d AC x r y r + dec yl) + — z r . (4.7) 

q q 

In a similar way as it was done for two-leg walks in the 
case of GM fractals (Appendix B) it can be deduced that 
do/q = b(b+ l)/2. Then, equation (|4.7| implies that the 
large r behavior of z r is governed by the larger of num- 
bers A 2 and b(b + l)/2. By computer enumeration of all 
possible one-leg configurations we found that inequality 
A 2 > b(b+l)/2 is satisfied for 2 < b < 8, so that z r - A 2r . 

With the established large r behavior of x r , y r and z r , 
and consequent behavior of A^ r \ and D^ r \ together 
with the previously found relation (|4. 3|) for B^ r \ it is not 
difficult to find out that scaling form for overall number 

(r) 

Z Q on MSG fractal is the same as in the case of corre- 



sponding GM fractal (|3 . 8[> . but with different values of 
uj and A. Consequently, the only exception - equality of 
the values of 7 for 3-simplex and SG fractals (6 = 2 case 
of MSG and GM fractals, respectively), seems to be ac- 
cidental. Indeed, this exception is also an indication that 
number of contacts per monomer is not the only relevant 
physical measure of HWs frustration. 

The particular values of the connectivity constant 10 
and the exponent 7 for MSG fractals with 2 < b < 8 
are given in Table HI whereas in Fig. [7] we depict 7 as a 
function of 1/6, for both GM and MSG fractals. One can 
see that 7 is not a monotonic function of 6, which has not 
been found for any of the critical exponents connected 
with SAWs on these fractal families HE H3, El- For 



n=3 



n=4 




0.0 0.1 0.2 0.3 0.4 0.5 

1/b 



FIG. 7: Values of the exponent 7 appearing in the scaling 



form Zo 



N~* found for the number of open HWs on 



GM (full triangles) and MSG (A) fractals, together with the 
corresponding values for SAWs (o) on the same fractal families 
[Tol | . as functions of the reciprocal of the fractal parameter b. 



the sake of comparison on the same figure we give the 
corresponding values of 7 for open SAWs [lfj , which are 
the same for GM and MSG fractals with equal b. In 
addition, one can also observe that values of 7 for HWs 
are separated in two groups: for even values of b exponent 
7 decreases when b grows, whereas for b odd 7 increases 
with b. This, however, resembles the behavior of the 
exponent ip, appearing in the scaling form X N fi N for the 
average number of different configurations of a branched 
polymer of N bonds on GM fractal with b > 2. Namely, 
for even and odd values of b in a recent paper [32j Dhar 
obtained two different exact expressions for ip. Finally, 
it is clearly seen in Fig. [7] that the difference between 7 
values for GM and MSG fractals becomes smaller with 
larger 6, so it seems plausible to investigate HWs only 
on MSG fractals (which is simpler than GM case), in 
order to obtain the large b behavior of the exponent 7. 
This limit is interesting because for large b already at 
the first stage of the construction of the GM fractal, one 
gets finite, but large, homogeneous triangular lattice, for 
which asymptotic form of the number of HWs has not 
been established yet. 



V. N-SIMPLEX LATTICES 

In order to extend our analysis to higher dimensional 
lattices, in this section we turn to n-simplex fractals with 
n > 3, which are embedded in d = n — 1 dimensional Eu- 
clidean spaces. To obtain an n-simplex lattice Q one 
starts with a complete graph of n points and replaces 




A, B, 



n=5 



B, A 2 B 2 A, 

n=6 




n=l 



B, 



B, 



-+ 



FIG. 8: All types of HWs needed for determining the connec- 
tivity constants and asymptotic forms of the numbers of HWs 
for n-simplex lattices with 3 < n < 7. Gray polygons denote 
n-simplex of order 



each of these points by a new complete graph of n points. 
The subsequent stages are constructed self-similarly, by 
repeating this procedure. After I such iterations one ob- 
tains an n-simplex of order I, which consists of Ni = n l 
points. The complete n-simplex lattice is obtained in the 
limit / — > 00. Fractal dimension df of n-simplex lattice 
is equal to df — Inn/ In 2. 



Closed HWs 

Any closed HW on n-simplex of order I + 1 can be 
decomposed into n parts within its n simplices of order I. 
Parts (steps) of the walk within the simplices of order I 
can be of [n/2] different types (see Fig. [5J, and we shall 
denote the numbers of corresponding HWs by B\, B2, 



Bu 



to 



72] 



The overall number of closed HWs is equal 



Z. 



E 



iiH K 



n 



B 



(5.1) 



where Ci 1 j 2 ... j., _ 1 .. 2] is the number of closed HW config- 
urations with ii, i[( n _i)/2] steps of B ly B[(„-i)/ 2 ] 
type, respectively. Numbers Bi satisfy a closed set of 
recursion relations: 



B 



0+1) 



E 



il «2-"t[n/2] 



ilH H[„/2]=™ 

l,2,-» J [n/2], 



[n/2] 

i=i 



(5.2) 



where b- '■ „■ is the number of HW configurations 
of .Bj-type that traverse the simplex of order (I + 1), and 
have i\, i[ n /2] steps of B±, B\ n /2] type, respectively. 
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Due to the self-similarity of the underlying lattice both 

(i) 

c iii2-i l (.n-i)/2] and Vi2 - i[„/2] depend only on n. Starting 



with the initial values B^ = 1, B, 



>(0) 



B 



(0) 
n/2] 



0, one can iterate recursion relations (|5.2p and calculate 

By' and in principle for any I, but since all these 
numbers grow very fast with I, it is useful to introduce 
new parameters: Xi = Bi/B[ n /2], i = 1, [n/2] — 1. 
Dividing relations for B^ +1 ^ K " RC+ 1 ) 
recursion relations: 



by B,,t' one obtains new 



(1+1) _ fi t^' 



,(») 

>/2]-l 



(0 

[n/2]-l 



(5.3) 



where 



fi ( x i ? 



l /2] 



- 1 )=^ +1) /« ) /2] ) 1 



-J - ^[n/2] / \ D 



i/2 



;(') 



V2] 



(5.4) 



i = 1, . . . , [n/2] — 1. Iterating these relations (starting 
with I = 1) one notices that large I behavior of x"p 
strongly depends on the parity of n. It turns out that 
for odd n values all xf tend to some finite constants, 
whereas for even n the following asymptotic relations are 
valid: 



which follows directly from (|5.4[) . Furthermore, from 
(15.71) one has 



InB 



(p) 



lnB<'> ^ 1 



(m) 



which, with p — > oo, gives 



In a; 



InB 



(0 



Y-L-\nf(x { { 



(m) 



x fc-lJ • 



m—l 



Then, since f{xi\ . . . sji-i) decreases with m, it follows 



„( m ) 



that n'lntj — const < In Ik. < rt'lnw, 
n lnw. Using this, together with (|5.6[) . for the numbers 
of closed HWs one finally obtains the scaling formulas 

f uj Ni , for n odd, 



InB 



(0 



z. 



(I) 

c 



for n even, 



(5.8) 



where a 



In 2 



t — — -T-, with df being the fractal dimension 

Inn df 1 J ° 

of the n-simplex lattice. Scaling forms and particular 
values for ui were found in [26] for n = 4, in [2^ | for 
?i = 5, and 6, and for n — 7 in Appendix [Cl of this paper. 



Open HWs 



,(0 



(v-y 



i = 1 , • • • , k — 1 , 



(5.5) 



where k — [n/2] and As is some finite constant less than 
1. Such behavior is a direct consequence of the fact that 
for even values of n none of Si, _B 2 , • • ■, -B/c-i-type con- 
figurations within a simplex of order (I + 1) can be ac- 
complished by n steps of B^-type through comprising n 
simplices of order I (see Appendix [Cj , which is not the 
case for simplices with odd n, because -B^-type steps oc- 
cupy even number of corner vertices of any simplex. 

With the established large I behavior of numbers Xi it 
can be shown that from (|5.1[) follows 



AW) 




for n odd, 
for n even. 



(5.6) 



Although these two asymptotic forms are not the same, 
from the definition (|1.2[) of the connectivity constant uj 
in both cases one gets 



muj = lim 



InB 



(0 



I— >oo n 



i+i 



lim ■ 

/ — too 



i/2] 



The last limiting value can be obtained via numerical 
iteration of the recursion relation 



InB 



(J+i) 



^Tr- = ^ + ^T^/(4°,-4'l 1 ) J (5-7) 



Any open HW on n-simplex of order (I + 1) can be 
decomposed into n parts within its n simplices of order 
I. The two ends of the HW are either both in the same 
or in two different Zth order n-simplices. In the first case 
the configuration of the corresponding HW part is of the 
two-leg type, whereas in the latter case the configura- 
tions of both ending parts are of the one-leg type. For 
each n there are n — [n/2] different types of one-leg con- 
figurations (for 3 < n < 7 they are depicted in Fig. [8J. 
The remaining (n — 1) or (n — 2) parts of the open HW 
correspond to B^-type steps. If we denote the numbers 
of one- and two-leg walks by and D^p , respectively, 
then the overall number of open HWs on n-simplex of 
order (I + 1) is equal to 



E^AfAf+^Df, (5.9) 

i,j=l 2—1 



where 



E, 



■ ' / j inV"i[ n /2] 

n/2] r— 1 

i r =n-2 



[n/2] 



n (*?>) ' 



[n/2] 



j ^ e n (^)"'( 5 - 10 ) 

[n/2] r=l 

J2 v=«-l 
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with coefficients e and / being the numbers of corre- 
sponding HW configurations. Numbers Af' of one- leg 
walks satisfy recursion relations which can be put in the 
following matrix form 



which, with obtained behavior of yp and £[„/ 2 ]-i' 
rectly gives: 



A 



a+i) 



A 



/ an 



V a p i 




(5.11) 



where p = n — [n/2] , and 



E 



x/2] 



[n/2] 

n 



B (i: 



(5.12) 



[n/2] 
r=l 



with a 



(y) 



being constant coefficients. Introducing 



»i »a—*[n/a] 

parameters ?/i = Ai/_B[ n / 2 ] ) * = 1) ■■■iPj an d dividing re- 
lations for A? +1) by B<£+J one gets recursion relations 



for the new parameters: 

/ y[ l+1) \ t 



\ 

where 



(i+i) 
Vp 



\ m p i 



Hp 



VP 



) 









(5.13) 



ay (Si, . . . , -B[„/2]) 



_ 9ij(xi,. . . ,x fc _i) 

" HJ /(xi, . . . ,x k -i)(B [n/2] ) n - 1 f(xi,...,x k -i) 

with <7y being polynomials of order n — 1, whose coef- 
ficients depend only on n, and / defined in (|5.4[) . Nu- 
merically iterating relations (|5.3p for Xi , and substituting 
them into functions m<j , after large number of iterations 
one obtains constant values, such that 



• for odd n, matrix m,j has only one eigenvalue A 
larger than 1, implying that yf^ 

• for even values of n, only m pp 



X 1 , 

1, whereas all 



other 



0, implying that y p tends to finite 

(0 



constant value, whereas all the other y\ tend to 0. 

Knowing the behavior of parameters y!p , one can get the 

large / behavior of the numbers A\ , and, consequently, 
the asymptotic form of the first sum in the right-hand 
side of (|5.9p . In addition, it can be shown that numbers 
of two-leg configurations Di are not necessary for estab- 
lishing the asymptotical behavior of (|5.9[) (terms with 
either scale as terms with A^Aj or they are much 
smaller). For n-simplex lattices with odd n all terms with 
A^Aj have the same asymptotical behavior, whereas 



for even n behavior of Z, 



(l+l) 



is governed by the term 

with (A p lJ ) 2 . In both cases the following asymptotical 
formula is valid: 



(i) 

'[n/2)-l 



D 



(I) 



i/2] 



di- 



; N ' N7 , 7 = 2^ , for odd 

/ ' ' Inn ' 



for even n. 



(5.14) 



Details of the derivation of these formulas can be found 
in Appendix [Cl Value of 7 = 2.5237... for n = 3 was 
found in Sec. II, whereas for n = 5 and 7, exponent 7 is 
equal to 2.1668... and 2.1079..., respectively. 

As in GM and MSG fractal cases, one observes that 
behavior of the number of long HWs on n-simplex lat- 
tices is strongly affected by the parity of the fractal pa- 
rameter. For lattices with odd n scaling forms are the 
same as for GM and MSG fractals. This scaling form 
coincides with the scaling forms for average numbers of 
TV-bonded SAWs on these lattices: with commonly ac- 
cepted symbols for critical exponents, average number of 
closed SAWs scales as [i N N a ~ 3 , whereas for the average 
number of open SAWs formula ^^iV 7-1 is valid (these 
are the same formulas as for SAWs on homogeneous lat- 
tices). Comparing these formulas with those obtained 
for HWs (and keeping in mind that formulas for HWs 
correspond to the overall numbers of HWs), one can say 
that corresponding value a — 2, obtained for closed HWs 
on all these lattices, certainly differs from the values ob- 
tained for SAWs BE! 03. The same holds for exponent 
7, meaning that HWs and SAWs belong to different uni- 
versality classes (see Fig. [7]). For HWs on n-simplex 
lattices with even n we obtained different scaling forms: 
Lo N {\ n J 2 ) N ° for closed, and ^(A^ 2 ^ 2 ')^ for open 
walks. However, one should note here that, whereas poly- 
mers on GM, MSG, and 5-simplex fractals for T > can 
exist only in swollen phase, on 4- (22j, and 6-simplex [24[ 
lattices, bellow some finite temperature Tg > polymers 
collapse into a more compact phase. Scaling of collapsed 
polymers has not been analyzed so far, and behavior of 
HWs might be an indication that on n-simplex lattices 
with even n, for T < Tg polymers do not scale in the 
same manner as they do for T > Tg. 



VI. SUMMARY AND CONCLUSION 

In this paper we have analyzed asymptotic behavior 
of the numbers of open and closed Hamiltonian walks 
on Given-Mandelbrot, modified Sierpinski gasket and n- 
simplex fractal families. Obtained scaling forms are sum- 
marized in Table HT1 One can see that scaling form (|1.1[) : 
Zn ~ u> N fig N a , proposed and obtained for some homo- 
geneous lattices with large number of sites N, is satisfied 
for all fractal lattices under consideration, but, with at 
least one of the exponents a and a being equal to 0. For 
HWs on GM, MSG, and n-simplcx lattices with odd n 
scaling form is ui N for closed, and oj n N' 1 (7 > 0) for 
open walks, i.e. exponent a formally may be taken as 
0. This value cannot be obtained with a simple gen- 
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TABLE II: Asymptotic behavior of the numbers of open and 
closed numbers of HWs on lattices with N sites. Values of oj, 
Xb, and 7 are given in the text, and d/ is the fractal dimension 
of the corresponding lattice. 



Lattice 


Closed HWs 


Open HWs 


all GM and MSG fractals 






n-simplex with odd n 






n-simplex with even n 


^ (ai) 





eralization of the formula a = (d — l)/d, valid for d- 
dimensional homogeneous lattice, but it can be explained 
following the original argument of Owczarek et al pj, 
which led to Ijl.ip for collapsed polymers on homoge- 
neous lattices. Namely, collapsed polymer chain on a 
homogeneous lattice has a form of a compact globule, 
with a sharp boundary separating it from the surround- 
ing solvent. Monomers on the boundary have smaller 
number of contacts with other monomers then those in 
the bulk of the globule, so that a factor /J,g s , with N$ 
being the number of monomers on the boundary, should 
appear in the scaling form. Furthermore, if one assumes 
that the boundary itself is homogeneous surface, then 
Ns ~ N^ d ~ x ^/ d , where N is the number of monomers 
in the polymer chain. The crucial point of this argu- 
ment is existence of two sets of monomers with different 
number of contacts. For HWs on both 3-simplex lattice 
and Sierpinski gasket (MSG and GM with 6 = 2, respec- 
tively), monomers at all sites, except for the three outer 
vertices of the lattices, have the same number of contacts 
(see Figs [T]and[3]), so that 'surface' factor with nontriv- 
ial a cannot arise. For GM fractals with b > 2 it can 
be shown that on the rth stage of fractal construction 
there are A^ r) = ^{[b(b + l)/2] r - 1} 'bulk' monomers 
with four contacts, whereas Ng — 3 = A^ r ) — Ng — 3 
'surface' monomers, with given by (|3.1[) , have two 

(r) 

contacts. This means that for r > 1 both N B and 

N ( g } scale as iV (r) ~ [b(b + l)/2] r , so that /if s can be 
incorporated into the term u> N , or formally one can put 
a = 0. For HWs on MSG fractals with b > 2, as well as 
on n-simplex lattices with odd n, all sites have the same 
coordination number (except for the finite number of cor- 
ner vertices of the whole lattice), and consequently all 
monomers have the same number of contacts. Hence, the 
argument of Owczarek et al for these lattices also gives 
(7 = 0. However, presence of the stretched-exponential 
term (A^ 2 )^"' for closed, and (a^ 2 )/ 2 ))^/ for 
open Hamiltonian walks on n-simplex lattices with even 
n, can't be a consequence of surface effects, since all the 
sites of these lattices also have the same coordination 
number. This means that there is some additional effect 
that should be considered in order to explain obtained 



scaling forms. It is not clear what that effect is, but one 
should notice that for Zth order n-simplex A/, 1 = 2', 
which is equal to the linear size of the lattice. Also, some 
clue about it can be achieved by careful inspection of the 
differences between the asymptotic behavior of the requi- 
site numbers of HWs for lattices with even and odd values 
of n. For instance, for odd n the numbers of long Bi-type 

HWs fulfil relation: B[ l) ~ B^ l) B?\ 2] , whereas 

for even n the number B^J 2 is much larger than any other 

Bf\ Hence, one can expect stretched-exponential terms 
in scaling forms for HWs on lattices on which 'entangled' 
configurations are more probable. This bears some anal- 
ogy with square lattice, for which it was found that large 
fraction of monomers participate in secondary structures 

To conclude, we can say that exact recursive technique, 
used in this paper for enumeration of all possible Hamil- 
tonian walks on Given-Mandelbrot, modified Sierpinski 
gasket and n-simplex fractal lattices, proved to be very 
efficient for open, as well as for closed walks. Recur- 
sion relations for the number of different types of com- 
pact configurations in simpler cases could be obtained by 
direct enumeration of the corresponding walks, whereas 
for lattices with more complicated structure we used ap- 
propriate computer programs. Once recursion relations 
were established, it was possible to generate and analyze 
numbers of very long HWs, eventually obtaining their 
asymptotical behavior. We believe that obtained scaling 
forms can help in better understanding of real compact 
polymers, therefore we intend to apply this technique to 
HWs on fractals embedded in three-dimensional space, 
in particular on 3d generalizations of GM and MSG frac- 
tals. It would also be interesting to extend the method to 
interacting HWs. Finally, we hope that our results might 
be useful in finding scaling forms for collapsed phase of 
interacting self-avoiding walks on fractals for which tran- 
sition from swollen to collapsed phase occurs at Tq > 0. 
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APPENDIX A: RECURSION RELATIONS FOR 
ONE-LEG HWs ON GM FRACTALS 

In this Appendix we prove that numbers of one-leg 
HWs: A, Ai, A2, and C, on Given-Mandelbrot fractals, 
satisfy recursion relations of the form: 
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A' = aAB m B[ l + a 1 A 1 B m+1 B[ 1 - 1 + a 2 A 2 B m+2 B^ 2 + cB m+1 B^ x C , 

A\ = (« ~ |) AB" 1 ^ 1 B™ +1 + a\A\B m B™ + (a 2 + p)A 2 B m+l Sf -1 + cB m B?C , (Al) 

A 2 = ^— ^A-B m ~ 2 B" +2 + — Aii? m_1 i?" +1 + (— + p) AsB" 1 ^ 1 - 1 + ~B m - x B? +1 C , 
2 2 V 2 / 2 



where p is the number of B-type (or -Bi-type) configura- 
tions within the (r + l)th order generator (see (|3.2p ). 
whereas a, 01, 02, and c are some integers depend- 
ing only on the scaling parameter b of the fractal, and 
m = b(b — l)/2, n = b — 1. In order to do that, we first 
notice that all configurations of -B-type can be divided in 
four sets, regarding the last step ("step" is here a part of 
the walk within the rth-order generator), as is sketched in 
Fig- EI a ) ■ The number of configurations within each set 
is denoted by rij, so that 711+712+713+714 = p. Each of the 
configurations from any of these four sets can be trans- 
formed into one A-type configuration by cutting the last 
step. In such a way, if the last step was -B-type it is con- 
verted into an A-type step, whereas a Bi-type step is con- 
verted into an Ai-type step. Number of configurations 
obtained in this manner, contributing to A^ r+1 \ should 
be doubled, due to symmetry (there are two corner rth 
order generators in which A- type walk can terminate). 
This means that (n 2 + 7i 4 )(BM) m+1 (B< r) )™ HWs of B- 
type on (r + l)th order generator can be transformed into 

2(7i2+?i4)(B w ) m ( J Bi r) )'M (r) HWs of A-type on the same 

generator. The remaining (ni + 7i3)(B^ r )) m+1 (B^ r '') ra 
HWs of -B-type on the (r + l)th order generator can be 

converted into 2(n x + n 3 )(B^) m+1 (B[ r) ) n - 1 A^ ) HWs 
of A-type (Fig. Ha)). 

In addition, an A-type configuration can have its last 
step in an interior rth order generator, and, as can be 
seen in Fig. [D(b) , that step can be of any of the four possi- 
ble one-leg types. Number of such walks with the A-type 
end is of the form n A {B^) k {B^) l A^ r \ where k and I are 
numbers of B and Bi-type steps, respectively. Numbers 
k and I must satisfy equation k + I + 1 = b(b + l)/2 — 2, 
because every rth order triangle must be traversed by 
HW, and also k + 21 + 2 = (b + 1)(6 + 2)/2 - 2, because 
every corner vertex of every rth order generator within 
the (r + l)th order generator, with the exception of two 
outer vertices, must be visited. The only solution of these 
two equations is k = m, I = n. In a similar way one can 
obtain the numbers of walks with the A\ , A 2 or C-type 
end in interior rth order generators, so that the number 
A(r+i) j s Gqua i to 

A (r+ ^ = (2n 1 + 2n 3 + n A1 )B m+1 Bl l ~ 1 A 1 
+ (2ra 2 + 2ti 4 + n A )B m B™A 
+ n A2 B m+2 B™-' 2 A 2 + n c B m+1 B'l- x C . 

On the other hand, from Fig. [^a) it is quite obvious that 
m = n 2 and 713 = 714, implying that n\ + 713 — 712 + 714 = 



I 

p/2, and consequently: 

A^+V = (p + n A )B m B[ l A+(p + n A1 )B m+1 B^- 1 A 1 
+ n A2 B m+2 B?- 2 A 2 + n c B m+1 B™- 1 C . (A2) 

In a similar way, both Ai and ^2-type configura- 
tions within the (r + l)th order generator can be di- 
vided in two classes: (1) configurations terminating in 
some interior rth order generator, and (2) configurations 
terminating in one of the two possible corner rth or- 
der generators. The number of Ai-type configurations 
in the first class is the same as in the corresponding 
case of A-type configurations, and the only difference 
between such A and Ai-type configurations is that B- 
step through one corner triangle should be replaced by 
Bi-step. These configurations then give rise to terms 
nj^B^B^A + n A1 B m B\A Y + n A2 B m+l B 1 ^ 1 A 2 + 

ncB m B"C in the relation for A^ +1 ' . The number of 
^2-type configurations in the first class can be obtained 
from the corresponding A-type configurations by dividing 
their number by two (because it is predefined in which or- 
der the corner vertices are visited) and substituting the 
two corner B-steps by B\ steps, thus leading to terms 
(n A B m - 2 B? +2 A +n A1 B m - 1 B? +1 A 1 +n A2 B m B™~ 1 A 2 

+n c B m - 1 B™ +1 C)/2 in the relation for A { 2 r+1) . 

As for the second class, lets first consider Ai-type 
configurations which visit the upper corner vertex (Fig. 
fTUT a)). and terminate in the right corner triangle. Each 
of them can be obtained from exactly one Bi-type con- 
figuration, by converting the last B or Bi-step through 
the right corner triangle into A or Ai-step, in a simi- 
lar way as was done in the case of A-type configurations 
(Bi-type configurations can be partitioned in the same 
way as B-type configurations in Fig. 0(a), with the only 
difference that B-step through the upper corner triangle 
is substituted by Bi-step). In addition, Ai-type configu- 
rations that visit the upper vertex can terminate in the 
upper triangle, implying that the last step is of A\ or 
A 2 type, whereas the step through the right-corner tri- 
angle is of type B. In Fig. [TUTb') it is shown how these 
configurations can be obtained from B-type configura- 
tions. Finally, one obtains recursion relation for Ai-type 
configurations in the following form: 

A[ = {^+n A )B m - 1 Bl +l A+{p + n A1 )B m B^A 1 

+ (p + n A2 )B m+1 B^- 1 A 2 + n c B m B[ l C. (A3) 

Since in the case of A2-type configurations the order 
in which corner vertices are visited is fixed, correspond- 
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(b) 

FIG. 9: (a) Possible configurations of B-type HWs on (r + l)th order generator of GM fractal, and corresponding A-type 
HWs, together with the corresponding numbers of walks. Small upward-oriented triangles represent the rth order generators. 
Encircled nodes are the corner vertices which are not visited by the corresponding HW. (b) Possible types of A-type HWs 
which terminate in interior rth order generators. 





(b) 



FIG. 10: Possible configurations of Ai-type HWs on (r + l)th order generator of GM fractal, which terminate in (a) the 
right-corner and (b) the upper corner rth order generator, together with the corresponding numbers of walks. 
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TABLE III: Coefficients a, eti, 0,2, c, and p, appearing in 
recursion relations (| Al|l and (|3.2p for the numbers of A, Ai, 
A 2 , B, and Si -type HWs on GM fractals with 2 < b < 7, 
found by computer enumeration of the corresponding HW 
configurations. 



b 


a 


ai 


a 2 


c 


P 


2 


2 


2 





2 


2 


3 


12 


14 


4 


22 


8 


4 


122 


128 


36 


212 


40 


5 


1842 


1532 


436 


2 704 


360 


6 


35 390 


23 812 


5 932 


42 368 


3 872 


7 


880 646 


486 284 


110 876 


878 168 


62 848 



TABLE IV: Coefficients a', a\, a' 2 , and c', appearing in re- 
cursion relation (|A5|I for the number of C-type HWs on GM 
fractals with 2 < b < 7. 



b 


a' 


oi 


a 2 


c 


2 


1 


2 


2 


6 


3 


12 


16 


8 


32 


4 


152 


168 


96 


352 


5 


2 544 


2120 


848 


4 048 


6 


52 072 


35152 


13 136 


67680 


7 


1 340 536 


735 312 


224 176 


1 374 944 



APPENDIX B: RECURSION RELATIONS FOR 
TWO-LEG HWs ON GM FRACTALS 

Here we prove that the numbers of two- leg HWs, D 
and Di, on GM fractals, satisfy recursion relation (|3.5[) . 
where coefficients fey are equal to: 

hi = p(m + l)B m B?, hi = 2pnB m+1 B^\ 
b 2 i = ^B m - 1 B? +1 ,b 22 =p{n+l)B m B?. (Bl) 

The two strands which form any two-leg configuration 
can terminate either in two different rth order generators 
within the (r + l)th order generator, or in the same one. 
In the first case, the corresponding number of walks is 
of the form B k B[XY, where X and Y are A, A\, A 2 , 
or C, and k + I + 2 = b(b + l)/2, since all the rth order 
generators must be visited. In the letter case the number 
of walks is of the form B k B[X, where X is D or D\, 
and k + I + 1 = b(b + l)/2. Each such configuration can 
be obtained by cutting some B or Bi-step of the B or 
Bi configuration within the (r + l)th generator. All the 
possibilities in which this can be done are sketched in 
Fig. [H so that Q5T) directly follows. 



APPENDIX C: RECURSION RELATIONS FOR 
HWs ON ^-SIMPLEX LATTICES 



ing walks from the class (2) can terminate only in one of 
the two corner triangles, lets say in the right corner one. 
Then, all possible cases, obtained from the partitioning 
of Bx-type configurations, are depicted in Fig. [TTJ These 
cases, together with configurations that terminate in in- 
terior triangles give the relation 



A > = , ^B m - 2 Bl l+2 A + ^^B m - 1 B\ i + 1 A 1 
2 2 1 2 1 

+ (— + p) B m B[ l A 2 + ^-B m - x B'l +1 C . (A4) 



Finally, with a = p + n^. ai = p + Uai, a 2 = n^, and 
c = nc from (|A2[) . (|A3[) , and (|A4[) . one gets recursion 
relations (|A1|) . Values of the coefficients a, a%, a 2 , and 
c for 2 < b < 7, are given in Table [TTT1 For 6 = 2 
and 3 we could find them by direct enumeration of the 
corresponding HW configurations, whereas for larger b 
we had to use computer facilities. 

In a quite similar manner it can be shown that the 
number of C-type HWs fulfills the recursion relation 



C = B m B1{a'^A + a[Ai + a' 2 ^-A 2 + c'C), (A5) 
B B\ 



where coefficients a', a' l7 a' 2 , and c' depend only on &, and 
their values for 2 < b < 7 can be seen in Table HV1 



In this Appendix we give recursion relations for the 
requisite numbers of HWs, needed for obtaining the scal- 
ing forms for the overall numbers of HWs on n-simplex 
lattices, for n = 4,5,6, and 7, together with some rele- 
vant details of the derivation of formulas (|5.8[) and (|5.14[) . 



4-simplex 

Closed HWs on 4-simplex were analyzed in [26j . Here 
we quote equations and results relevant for our present 
analysis. Recursion relations, as well as the initial values 



of the numbers B± and B^ are 



B' 



2B\ + AB\B 2 + 6BfB 2 2 , b[ 1] 

(i) _ 



B' 2 = B$ + AB\B 2 + 22BI , B [ 2 L> = 1 



(CI) 



The requisite number x = x\ = B\fB 2 satisfies recursion 
relation 



B[ _ 2x 4 + Ax 3 + 6x 2 
~W 2 ~ x A + 4x 3 + 22 



2. 



(C2) 



Explicit numerical iteration shows that x^ — > 0, when 
I — > oo, so that recursion relation obtains the approxi- 
mate form 



,(*+!) 
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.(0 
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2n l B'B"A 2 nS'-'B'C'A, 2n 3 B"B"A 2 n^B^A, 

FIG. 11: Possible configurations of Tl^-type HWs on (r+ l)th order generator of GM fractal, with the ending step in the corner 
rth order generator. 




pBT%" p B" +, b: pB"B"' pB"B l 




(m+l)pB"B"D 2npB l " 1 B;'-'D l mpB ml B 1 n+i D 2(n+l)pB m B 1 n D 1 

FIG. 12: All possibilities of cutting B or Bi-type configurations within the (r + l)th order generator of GM fractal, leading to 
all possible D or Di-type configurations, in which the two legs terminate in the same rth order generator. Terms under each 
configuration are equal to the corresponding walks, within the (r + l)th order generator. One should notice that by cutting 
steps of B\ configurations, due to symmetry, doubled number of Di configurations is obtained (since it is predefined which of 
the two legs visits the upper corner vertex). 



This relation implies that x^ ~ A^ , and precise numeri- 
cal analysis gives As = 0.836620 . . .. The overall number 
of closed HWs on 4-simplex of order (I + 1) is equal to 



= 3[B\ 



>(0 



which, with established asymptotic behavior of x"\ for 
large I gives the same form as in (|5 . 6|) : 



1} - (a« } ) 4 . 



and, consequently, Zq ~ uj Nl (\g) N i , with lo = 
1.39911.... 

Directly enumerating one-leg HW configurations, we 
have found that numbers A\ and A 2 satisfy recursion 
relations of the form (|5.1ip with 

a n = 6B?0Bi + B 2 ) , a 21 = B 2 1 {2B l + 3B 2 ) , 
a 22 = 22B 2 (B 1 +B 2 ) + B 2 (7B 1 + 16B 2 ), 



and a± 2 = 60,21, so that for parameters y\ = A\j B 2 , 
y 2 = A 2 /B 2 , one obtains relation (|5.13[) . where 

6x 2 (1 + x) x 2 (3 + 2x) 



and 



m 22 = 



22 
22 



4a: 3 - 
22a: 



x^ 

■ 16a; 2 



m 2 i 



lx A 



22 + Ax 3 



22 + 4x 3 + x A 

?71i2 = 6m21 



Since x^ — B^ / B 2 tends to 0, it is obvious that mn, 
777,12, and mai tend to 0, whereas 77122 —> 1, when I — » 
00, which implies that y\ tends to 0, and y 2 to some 
finite constant value. Explicit simultaneous iteration of 
relations (|C2[) and (|5.13[) . starting with the initial values 

x^ = 2, — 6, and y^ — 2, indeed quickly shows 

that y{ l) -► 0, and yf -> 62.1081 . . .. 

The overall number of open HWs on (I + l)th order 
4-simplex is equal to 



= 12B 2 (A 2 + 2A 1 A 2 + 3A 2 + B 1 D 1 ), (C3) 

where superscript (I) was suppressed on the right-hand 
side, and D\ is the number of one of two possible two-leg 
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FIG. 13: Possible types D\ and D2 of two-leg HWs on 4- 
simplex lattice of order I. 

HWs (see Fig. [T3|) on the Zth order 4-simplex. To obtain 
asymptotical behavior of the numbers D\ and D 2 , one 
needs recursion relations: 




dn = \2B X B\ + \2B\B 2 + 8B 3 , 

d 21 = 2Bf(Bi + 3B 2 ), d 22 = AB\ + 88B*, 

du = 4d 2 i, and Fi = A 2 1 B 1 (<iB 2 +6B 1 )+24:A 1 A 2 B 1 (B 2 + 
Bi) + A\ B 1 {A2B 1 + 64B 2 ), F 2 = A\B\ + §A X A 2 B\ + 
Aj(l6B 2 + AABiB 2 + 66B|). Introducing new variables 
wi = D\/B 2 and w 2 — D 2 /B 2 one gets new recursion 

relations, which show that wf — > for I — > 00. Since 
from Eq. (IC3|l follows that 

= 12a; 2 B 4 (y 2 + 2yi y 2 + 3y 2 + x Wl ), 

one finds that two-leg HW configurations do not con- 
tribute significantly to the overall number of open HWs 
for large I, and, furthermore, that the following asymp- 
totic formula is valid: 

4 Z+1) -M 2 K) 4 - 

Finally, since B<p ~ io 4 \ and ~ A^, one obtains 

zW ~ , \ N ' /2 

Z U X B ■ 



TABLE V: Coefficients a™' given by (|5.12|l . which appear in 
the recursion relations (|5.11|) . for 5- and 6-simplex lattices. 



5-simplex lattice 



h 


«2 




11) 
a: - 
«i « 


(21) 

op .' 

•"2 


(22) 


(31) 

»1 


H21 

Ol 

11 »a 


(sal 
a l33) 





4 




24 


44 


548 


22 


428 


1042 


1 


3 







44 


636 





296 


472 


2 


2 




72 


38 


374 





104 





3 


1 




72 


24 


140 


4 


36 





4 










6 


26 


1 


6 


5 






Q 


(12) _ l 9fv (21) 
ilia _ Wa i,is' 


a (13) - 24a (31) 


(23) 

a = 


r, (32) 










6- 


simplex 


lattice 






ii 


«2 


«3 


(11) 

il »aia 


(21) 

a- ■ • 


(22) 
OL- ■ ■ 
1713 


(31) 
OL- ■ ■ 
ii 1713 


(32) 

a ■ ■ ■ 

"»1 W3 


Q (33) . 

"»1 W.*3 








5 

















541568 





1 


4 

















541568 





2 


3 








57600 





47168 


447232 





3 


2 








82848 





56768 


323264 





4 


1 


2640 


2568 


66744 


1042 


35256 


148000 





5 





2640 


1776 


23688 


614 


10220 


33160 


1 


2 


2 








25008 





9600 


47168 


1 


3 


1 





1776 


45552 


472 


17352 


58464 


1 


4 





2640 


2172 


27768 


538 


9576 


25586 


2 


1 


2 








4248 











2 


2 


1 





624 


14616 





3460 


6768 


2 


3 


n 
U 


2280 


1350 


15360 


192 


39 10 


77HO 

1 /U8 


3 





2 








456 











3 


1 


1 


480 


216 


2928 





480 





3 


2 





1800 


612 


5358 


68 


1060 


1188 


4 





1 


120 


36 


300 


5 


56 





4 


1 





720 


180 


1182 


23 


206 


168 


5 








120 


24 


126 


3 


20 


21 


Q 


(12) 
»1»2 


>3 


on (21) (13) , on (31) 


a (23) 
3' "n*2* 


_ c ( 32 ) 
3 — n^3 



5-simplex 



According to 27] recursion relations for requisite num- 
bers B\ and B 2 , and their initial values, needed for ob- 
taining the number of closed HWs are such that x^ 1 ' = 
B (i+i) = b { 1 +1) /B { 2 1+1) = /i(ar»)//(a:W), x 0) = 3, with 



fi 
f 



6a; D 
2x 5 



30a; 4 
13x 4 



78a; 3 
32a; 3 



96a; 2 + 132a; + 132, 
88a; 2 + 220a; + 186. (C4) 



Numerical iterations give — > 0.8023188..., whereas 
further analysis of relation (|5.7p reveals that 10 = 



1.717769.... and B 



Since 



7(1+1) 



12BI + 30B*B 2 + 60BlB% + 132B| . 



it is obvious that for large I one obtains Z^ 1 ^ 



B, 



(0 



and consequently Z { 



(i) 



For 5-simplex lattice there are three types of one- 
leg configurations, A±, A 2l and A3, which, as we have 
found by computer enumeration, fulfil recursion relations 
([5TTT]) . with coefficients af^ (see Eq. frig]) ) given in Ta- 

ble[Vj Parameters yf* = Ar* /B% satisfy matrix relation 
(|5.13p , with all coefficients my tending to finite constant 
values when I — ► 00, such that corresponding matrix has 
one eigenvalue larger than one: A = 5.718.... This means 
that yf^ ~ A', which was also confirmed by direct itera- 
tion of the recursion relations for x and yi. 

Overall number Zq +1 ^ of open HWs is of the form 



(|5.9p . and, similar to the 3-simplex case, it turns out 
that terms with the numbers Di of two-leg HWs have 
the same large I behavior as the terms with numbers 
of one-leg HWs. Consequently, explicit quoting of re- 
cursion relations for Di numbers is not necessary here, 
and for the number of open HWs one obtains: Zq 1 ^ ~ 
{B^ fA^Af - \ 2l (B% ] ) 5 . Then, straightway follows 
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that ~ oj n 'N?, with 7 



2 In A/ In 5 = 2.1668. 



6-simplex 



From the recursion relations for the numbers B± , 



B$\ and £ 3 . 



?(') 



obtained in [27 



certainly follows 

that xr x> = fi{x y ;>,x^)/f(x y ;>,x y 2>), and x { 2 +1) = 
h{x[ ,x 2 l) )/f(x\ l \x 2 with A, / 2 , and / being poly- 
nomials of order 6. Iterating these relations one con- 



cludes that both a 

M ^ J 1 ) 



0) 



and 



„0) 



1 , auu ^ 2 tend to 0, and also 
x\ l> x^ 1 , so that for large I the following approximate 
relations are valid: 



„0+i) 



25008 ( (Q,4 
541568 V 2 ; 



„0+D 



94336 
541568 



(0\2 



From the latter of these relations follows 



,(') 



X 2 



and then, from the first relation one gets x\ ' 

(A 2 ^) . Precise numerical analysis gives Xb 
0.9888.... Analysis of the corresponding relation ([57 

' ,6 ' Finally, since 



U) 



gives uj = 2.0550047..., and B 3 
3960B 2 B| + 792051^1 + 5580^, for large I one ob 

6 / r „,\6 



117QBfB'l + 19205^1 



tains Z, 



0+1) 
c 



const 

i.ln2/ln6 



and 



Z«~^ ! (A|)(«') 

Numbers A±\ A 2 , and v4 3 '' of one-leg HWs for 6- 
simplex lattice satisfy recursion relations (|5. 1 1|) with co- 
efficients af^) i given in Table IVl Direct numerical it- 
eration of the recursion relations for xi and yi parame- 
ters, quickly shows that and y 2 tend to 0, whereas 

y^ approaches large, but finite constant: 4. 08311. ..10 5 . 
This is in accord with the fact that 77133 is the only ele- 
ment of the matrix m of the relation (|5 . 13[) which for 
large I tends to finite constant, i.e. not to 0, as do 
all the other my. Finally, the overall number (|5.9p of 
open HWs is for large I governed by the term A 3 B 2 , 



so that Z, 

7(0 . 



O+i) 



(4°) (4°) - 



and 



(A|) 



(6 



IUn2/ In 6 



7-simplex 

To enumerate closed HWs on 7-simplex lattice one 
needs three requisite numbers B\, B 2 and B 3 (see Fig. 
[5]), which satisfy recursion relations of the form 



?«3 

^3 ' 



1,2,3. 



i\— 22— «3— 



Coefficients &^ 2 j 3 > found by explicit computer enumera- 
tion of the corresponding HW configurations, are given 



TABLE VI: Coefficients appearing in recursion relations B[ = 

J2 b i%2i3 Bi i B ^ B 3 3 for the numbers Bi (i = 1,2,3), needed 
for enumeration of closed HWs on 7-simplex lattice. 



h 


42 


13 


, (I) 


b-' ■ 

»1 »9«H 


1 (31 








7 


64988160 


94599168 


115468800 





1 


6 


64988160 


141735936 


165548544 





2 


5 


53667840 


118798848 


106301952 





3 


4 


50112000 


78990336 


49499520 





4 


3 


32536320 


38795520 


19747584 





5 


2 


14097600 


12925632 


5819328 





6 


1 


3787680 


2712192 


1077136 





7 





473760 


275616 


96864 


1 





6 





6498816 


11372928 


1 


1 


5 





10733568 


9748224 


1 


2 


4 


5660160 


15033600 


6708480 


1 


3 


3 


9319680 


13014528 


5011200 


1 


4 


2 


7985040 


7048800 


2440224 


1 


5 


1 


3459360 


2272608 


704880 


1 


6 





626400 


331632 


94692 


2 


1 


4 





566016 





2 


2 


3 


812160 


1397952 


283008 


2 


3 


2 


1953600 


1597008 


349488 


2 


4 


1 


1449360 


864840 


199626 


2 


5 





396720 


187920 


43242 


3 


1 


3 





54144 





3 


2 


2 


266880 


195360 


20304 


3 


3 


1 


384000 


193248 


32560 


3 


4 





163800 


66120 


12078 


4 


1 


2 


33600 


13344 





4 


2 


1 


80400 


28800 


3336 


4 


3 





50400 


16380 


2400 


5 





2 


3720 


672 





5 


1 


1 


12240 


3216 


336 


5 


2 





11520 


3024 


402 


6 





1 


840 


204 


31 


6 


1 





1680 


384 


51 


7 








120 


24 


3 



in Table IVH whereas initial values are: b[^ = 120 



(i) _ 



_> = 24, and B^ = 3. Introducing parameters 



B 

Xi = B1/B3, X2 = B2/B3, one obtains closed set of 
recursion relations 



_ , _ fl(xi,X 2 ) , _ f2(xi,X 2 ) 



f(xi,x 2 ) 



f(xi,x 2 ) 



(C5) 



where 



7 7 

i 1= i 2 =0 
7 7 



'2 

2 • 



i 1= i 2 =0 
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f{x 1 ,X 2 ) = E E 



i 1= i 2 =0 



These recursion relations quickly lead to the conclu- 



sion that 



0.690015. 



., and x^p 



1.14695. 



as i — > oo. The limiting values are indeed one of 
the solutions of the algebraic system of equations x\ = 
fi(x 1 ,x 2 )/f(xi,x 2 ), x 2 = f 2 {x\,x 2 )/ f{x\,x 2 ), actually 
the only positive one. 

The overall number of closed HWs on (I + l)th 

order 7-simplex is equal to 

7 7 



z c +1) -EE <vt , ,.n\n^ ■ j 



EE 



i 1= i 2 =0 

which for large I obtains approximate form: 



AW) 



const 



In 5 



(i+i) 



In B 



(i) 



+ ^ln/(^ ) ,4 Z) ), 



ji+i 



V 



(C6) 



which can be numerically iterated, together with the 
recursive relations |C5j for x± and x 2 . In that way, 



one finds In a; = limj^oo {\\\B^ /!' 



0.87382 . 



i.e. 



uj = 2.396056. 
write 



In addition, from relation (|C6[) one can 



InB 

Jk 



(fe) 



hxB. 



?+E^iW(4 ro) ,4 m) ): 



so that when k — » oo follows 



lnS 3 

mcj = 



(/) 



Since f(x\ 
is less then 



7' 



(m) (m) 



("0 2;( m ^ 



) decreases with m, the above sum 



<i) 



6*7' 

5.1, 



- , and consequently 7 l In w > 
for all I. This means that for 



I 1 one has m-B^ 



7' lnc 



B 



(0 



finally 



implying 



Recursion relations for the numbers of four pos- 
sible one-leg HWs are too cumbersome to be quoted 
here, but are available upon request. Numerical analysis 
of corresponding relations l|5.13p for parameters = 



A 



(l) /B (D 



reveals that y. 



(i) 



X\ with A = 7.7749..., be- 



meaning that InZp /7 l+1 ~ InB^ /7 l . From the re- m ^ ne limit I 
cursion relation for Z?g follows equation: 



ing the largest eigenvalue of the matrix rriij , obtained 



oo. All terms in the formula (15.911 for 



the number Zq +1 ' of open HWs are of the same order of 
magnitude for large I, so that 



and, consequently: 
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A') 



(7') 



In A 
! h77 



2.10791. 
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